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"Glassy Dynamics" in Ising Spin Glasses — Experiment and Simulation 
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The field-cooled magnetization (FCM) processes of Ising spin glasses under relatively small 
fields are investigated by experiment on Feo.55Mno.45Ti03 and by numerical simulation on the 
three-dimensional Edwards- Anderson model. Both results are explained in a unified manner by 
means of the droplet picture. In particular, the cusp-like behavior of the FCM is interpreted 
as evidence, not for an equilibrium phase transition under a finite magnetic field, but for a 
dynamical ('blocking') transition frequently observed in glassy systems. 
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Slow relaxational dynamics of noncxponential type has 
been observed in various disordered systems, and is now 
called 'glassy dynamics'.^ The most interesting slow dy- 
namics is the one observed in systems, such as structural 
glasses and spin glasses, whose relaxation processes are 
determined cooperatively through interactions between 
constituent elements of the system. In this case, the dis- 
tribution of relaxation times exhibits a peculiar temper- 
ature dependence besides the one attributed to the ther- 
mally activated process. A typical example is the crit- 
ical slowing down in Ising spin glasses, which are con- 
sidered to possess an equilibrium phase transition at a 
finite temperature, Tc, under zero magnetic field. On the 
other hand, a glass transition temperature, Tg, of a struc- 
tural glass is considered to be a 'blocking' temperature in 
the sense that relaxation modes, whose relaxation time 
significantly exceeds the experimental time window, ap- 
pears to be frozen below Tg. By the word 'glassy dy- 
namics', here, we mean such slow dynamics of a coop- 
erative origin which, at the same time, involves thermal 
'blocking' at low temperatures. We argue in the present 
letter that the field-cooled-magnetization (FCM) process 
in Ising spin glasses is one typical example of such 'glassy 
dynamics'. 

In spin-glass (SG) studies, one of most notable prob- 
lems yet unsettled concerns the stability of the equilib- 
rium SG phase under a finite magnetic field h? This is 
the case even for Ising spin glasses which are the con- 
ceptually most simple SG systems. The mean-field the- 
ory predicts the stability of the SG phase up to a cer- 
tain critical magnitude of h specified by the de Almaida- 
Thouless (AT) line,'^ while according to the droplet the- 
Qj,y^4-6 ^YiQ SG phase is unstable even in an infinitesimally 
small h. Recently, Takayama and Hukushima (TH)'^ have 
carried out numerical simulations of the field-shift aging 
protocol on the three-dimensional (3D) Ising Edwards- 
Anderson (EA) model, the results of which strongly sup- 
port the droplet picture. They have claimed that the AT- 
transition-like behavior observed experimentally^ is well 
interpreted as a dynamical crossover from the SG behav- 
ior to the paramagnetic one upon application of a field 
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h!^ Further numerical works in favor of the droplet pic- 
ture have been reported more recently.^"' 

Since early SG studies, it has been observed that the 
zero-FCM (ZFCM) deviates from the FCM at a certain 
temperature, denoted here as Xirr, and that the FCM ex- 
hibits a cusp-like behavior at a certain temperature, de- 
noted here as T* , below which the FCM becomes nearly 
independent of T. Here, the ZFCM (and FHM used be- 
low) is the magnetization measured in a reheating pro- 
cess after a ZFC (FC) process. Without consideration of 
cooling/heating rate effects, these phenomena have often 
been considered to be evidence for the SG phase tran- 
sition at Till- — T* under a finite field. The main pur- 
pose of the present study is to address this problem by 
performing detailed experiments on the Ising spin glass 
Feo.55Mno.45Ti03 with ~ 22.3 K,!^ by extending the 
Monte Carlo simulation on the same Gaussian Ising EA 
model as in the study of TH with ~ Q.^h^^^^^ (in 
units of the variance of interactions), and by comparing 
the results of the two analyses. The consequence is that 
the above-mentioned FCM behavior, a most fundamen- 
tal phenomenon in spin glasses, can in fact be interpreted 
as the 'glassy dynamics' in the sense described above. 

Let us begin with two sets of FCM, FHM and ZFCM 
curves, measured with two different cooling(heating) 
rates, and shown in Fig. 1(a). It is clear that T{„, the tem- 
perature at which the irreversibility between the FCM 
and ZFCM occurs, does depend on the cooling rate. Fur- 
thermore, is significantly higher than 
T*, the cusp-like temperature of the FCM.^^ The corre- 
sponding FCM and ZFCM obtained in our simulation are 
presented in Fig. 1(b). The rate-n in the figure means the 
cooling process with n MC steps (mcs) at each T, which 
is changed in steps of AT = 0.01. When the rates are 
measured in units of Tc/^o, with Iq being the microscopic 
spin flip time (which is 1 mcs for simulation and 10^^^ s 
for experiment), the rate of the process with rate-100 in 
the simulation is 1 • 10""*, which has to be compared, 
for example, with 2 • 10~^^ of that in the experimental 
process denoted as 1.7 mK/s. The results in Fig. 1 im- 
ply that the occurrence of the irreversibility between the 
FCM and ZFCM corresponds by no means to a certain 
equilibrium phase transition. 

At T > Tirr, the FCM, which coincides with the 
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Fig. 1. Cooling/heating rate dependences of FCM, FHM and 
ZFCM of (a) Feo.55Mno.45Ti03 with h = 5 Oe, and (b) EA 
model with h = 0.1 and rate-1, 3, 10, 33, ... 3333, and 10000 
(without FHM). The arrows in (a) indicate the cooHng/heating 
measurements. The Une in (b) represents the Curie law. In the 
inset of (b) the enlarged FCMs with rate-333, 1000, 3333 and 
10000 (from top) are shown. The averages were taken over at 
least 3800 samples for the data used in the present work. 
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Fig. 2. The FCMs (FHMs) observed in processes with a stop 
of the cooling at Tstp for a period fgtp are shown by thick 
solid (dashed) curves (ft = 5 Oe and 0.1). (a) The results for 
Feo.55Mno.45Ti03 with Tstp =22.8 K, 22 K and 19 K (top to 
bottom). Here, the upper and lower curves are shifted in height 
for clarity, and the FCMs (FHMs) without the halt of cooling 
are also shown by thin solid (dashed) curves. The inset of (b) 
shows the time evolution of FCM during the stop. 



ZFCM, is the magnetization M in equilibrium under 
ft- > 0. The small deviation of this equilibrium M/h from 
the Curie law may be attributed to the nonlinear field 
effect in the simulation with /i = 0.1 (see Fig. 4(b) be- 
low), and to a small nonvanishing Curie- Weiss constant 
for the experiment with h = 5 Oe. Here, we note that 
h = 0.1 in the simulation corresponds to the order of 1 T 
in the experiment on Feo.ssMno. 45X103. In this tempera- 
ture range, the SG short-range order, represented by the 
SG coherence length £,{T), and the corresponding corre- 
lation time, t{T), increase, and at T ~ Tii-r the latter 
reaches the time scale proportional to the inverse of the 
cooling rate. Actually, if this proportionality constant is 
set to unity, and Ti„ is estimated using the peak temper- 
ature of the ZFCM, our simulation results can be fitted 
to an expression of the critical slowing down 

r(r)^ao[(r-rg(ft))/rg(M]— , (i) 

with gq ~ 30, zi^ = 11.0 and Tg{h = 0.1) ~ 0.81, where 
Tg{h) is a supposed transition temperature under a fi- 
nite h. Here, we don't claim that Tg(h) is an equilibrium 
transition temperature. Instead, we want to emphasize 
that the SG short-range order starts to increase at tem- 
peratures as high as 2Tc, with the value of zi^ compara- 
ble to that obtained for the equilibrium transition under 



h — 0.^'^'^^ We therefore consider that the SG correlation 
of the order of ^* ^ T^^^{Ti„) has appeared at T ~ Tirr 
on average, and that locally correlated regions of rela- 
tively large size start to be thermally blocked around 
Tin-. We call them spin clusters. 

At T < Tirr, not only the ZFC state but also the FC 
state is out of equilibrium. This is clearly seen in FCM 
processes, where the cooling is stopped at a certain tem- 
perature, denoted as Tgtp. As shown in Fig. 2, the FCM 
increases during a stop made in the temperature range 
Thr > Tstp > Tc (Tstp = 22.8 K in the experiment and 
Tstp = 1.2 in the simulation). Because of a fairly large 
dynamical exponent z in eq. (1), the size of the SG short- 
range order is not expected to appreciably increase fur- 
ther during the stop within this temperature range. In- 
stead, spin clusters thermally blocked at around Tirr tend 
to equilibrate, or further polarize, yielding the FCM in- 
crease observed both in the experimental and simulation 
results. In an FC process without such a stop, spin clus- 
ters of a smaller size than the already blocked spin clus- 
ters are blocked, and at around T* , these blocked clusters 
are considered to come into contact with each other. 

The observed FCM behavior at temperatures compa- 
rable to and lower than T* can be interpreted in terms of 
the droplet picture,''"^ where the field crossover (corre- 
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lation) length,^ Lh, plays important roles. separates 
the characteristic behavior of droplet excitations in the 
equilibrium SG phase below Tc by their size L; it is dom- 
inated by the Zeeman energy (~ y/q^AhL^^^) for L > 
and by the SG stiffness energy (~ TL^) for L < Lh- Here, 
qea is the SG order parameter, d the spatial dimension, 
T(r) the stiffness constant of SG domain walls, and 6 
the gap exponent. Explicitly, Lh is written as 

U ^ {T{T)/h^f ^ [((T, - T)lT,r^"/h]\ (2) 

where 5 = (c?/2 — 0)~^ and the exponent acs (> 0)^' is 
given by the temperature dependence of T and (/ea- 

We now apply this droplet picture to spin clusters 
grown in the nonequilibrium FC process. At tempera- 
tures higher than r*(w Tc), spin clusters are thermally 
blocked independently of each other, or in other words, 
the SG stiffness energy introduced above is not effective 
at all. Thus, the configurations of the spin clusters are 
overbalanced to the Zeeman energy. At and below T*, 
the SG stiffness energy becomes effective and competes 
with the overbalanced Zeeman energy. We then natu- 
rally expect reorientation of spins, which results in the 
growth of new spin clusters with SG short-range order for 
a given set of (T, h) within old spin clusters.^* We denote 
the mean size of the new spin clusters simply as ^. The 
consequences of this restoration of the SG stiffness en- 
ergy, or the growth of i^, are not only the FCM cusp-like 
behavior at T* but also a decrease in FCM by further 
cooling the system below T* . For the same reason, an 
initial decrease in FCM is expected during a stop of the 
FCM process at Tgtp below T* . All of these expectations 
were in fact observed in the experiment as can be seen in 
Figs. 1(a) and 2(a). In the simulation, the initial decrease 
in FCM at the stop at Tgtp = 0.8 (< Tc) is clearly seen 
in Fig. 2(b), and the decrease in FCM below T* during 
cooling can be recognized in our slowest-cooling process 
with rate-10000 shown in the inset of Fig. 4(b) below. 
This FCM decrease below T* is contrary to the FCM 
increase commonly observed in noninteracting magnetic 
nanoparticle systems. "'^^ 

During a stop at a temperature below T*, the SG 
short-range order of size ^ continues to increase and so 
the FCM decreases until ^ becomes comparable with Lh 
given by Eq. (2). Then, ^ no longer increases since the 
system is now in the paramagnetic state. The magneti- 
zation makes a turn and increases (due to the polariza- 
tion of spin clusters) up to the equilibrium paramagnetic 
value. This is the dynamical crossover scenario proposed 
by TH. However, for ft, of a few Oe, whose corresponding 
Zeeman energy in Feo.55Mno.45Ti03 is about lO"* x Tc, 
Lh is in general much longer than and it will take 
an astronomical amount of time for ^ to reach Lh- An 
exception is expected at Tgtp close to Tc, at which Lh can 
be small (Eq. (2)) and the growth rate of ^ is the fastest. 
Actually, we experimentally observed the expected FCM 
upturn behavior at such a stop as shown in Fig. 3. Sim- 
ilar behavior has been observed for Feo.sMno.sTiOa,^^ 
while it is vaguely seen in the simulation with rate- 100 
at Tstp = 1.05 (see inset of Fig. 2(b)). 

Another finding observed in both Fig. 1(a) and 
Fig. 4(b) is that the FHM is smaller than the FCM. This 
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Fig. 3. Time evolution of FCM during a stop of the cooling for 
Feo.55Mno.45Ti03(/i = 5 0e). 
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Fig. 4. FCM and ZFCM curves under different fields and with 
different cooling rates for Feo.55Mno.45Ti03 (a), and EA model 
with rate- f 0000 (b). fn (b) the corresponding FCM curves with 
rate-333 are shown by thin lines and in the inset the FCM and 
FHM curves are shown at an enlarged scale. 



can be simply regarded as the cumulative memory effect 
of the SG short-range order. ^-'^ The SG order represented 
by ^ increases at T below T* and so the depolarization 
increases but ^ is still much shorter than Lh when the 
system is heated back to temperatures near T* . Also one 
can see that the FCM at lowest temperatures is smaller, 
the slower the cooling rate. In the simulation, this is ob- 
served in the cooling processes slower than rate-333, as 
indicated in the inset of Fig. 1 (b) . In the experimental re- 
sults shown in Fig. 1(a), crossing of the two FCM curves 
is observed. A possible interpretation for this is that, 
for the slower cooling process, ^* is longer and the cor- 
responding depolarization effect due to the SG stiffness 
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energy below T* is larger. 

Lastly, let us discuss the /i-dependences of the FCM 
and ZFCM curves. For the experimental results shown 
in Fig. 4(a) one can notice the following characteristics. 
As h increases from 5 Oe, which we have thus far inves- 
tigated, Tirr approaches T*. This implies that a stronger 
field makes it easier for the ZFCM curve to merge with 
the FCM curve in the paramagnetic phase. Under a suf- 
ficiently large h, on the other hand, Ti„ becomes lower 
than T* , whose associated cusp- like shape is much more 
rounded. The phenomena in this field range coincide with 
those observed for Feo.sMno.sTiOa by Aruga Katori and 
Ito,* and are interpreted as the dynamical crossover by 
TH; ^ in the ZFCM process reaches a relatively shorter 
L?i at T « Tirr within the observation time at each 
temperature. In the processes under /i of a few hun- 
dred Oe, we see Tjrr ~ T*. A rounded cusp-like shape 
in the FCM curve in these processes can be regarded as 
a blocking phenomenon since a difference between the 
FCM curves with different cooling rates is still notice- 
able below T* though it decreases with increasing h. As 
shown in Fig. 4(b), Tirr obtained by the simulation also 
decreases significantly as h increases. 

We have argued that the FCM behavior of Ising spin 
glasses observed in both the experiment and simulation 
can be interpreted in a unified way by a scenario based on 
the (extended) droplet picture, thereby we have empha- 
sized the viewpoint of 'blocking' of the FCM, or 'glassy 
dynamics'. The scenario proposed in the present work is 
rather intuitive and qualitative, but, at least at low tem- 
peratures under relatively large h, it is consistent with 
the results obtained in our related works,^'^^ where we 
have argued that the experimental and simulation re- 
sults, whose timescales in units of microscopic spin-flip 
time differ by more than ten orders of magnitude, can 
be understood in a unified way even semiquantitatively. 
Combined with our most recent work claiming the dy- 
namical breakdown of the SG state under relatively small 
h,^^ the present results strongly suggest that the slow 
dynamics in the FCM processes we have observed under 
moderate h is 'glassy dynamics' far from equilibrium. 

The viewpoint of 'glassy dynamics' is a natural con- 
sequence of the marginal stability of the SG phase pre- 
dicted by both the droplet theory and the mean-field 
theory. It implies an infinitely wide distribution of re- 
laxation times, at least mider h = 0. In the 1980's, the 
peculiar FCM phenomena discussed thus far had already 
been observed in various spin glasses and were reviewed 
by Lundgren et al.^^ They introduced a finite equilibra- 
tion time, teq, at all temperatures; however, the above- 
mentioned marginal stability implies an infinite t^q at 
T < Tc under h = 0. According to the droplet theory 
teq is expected to become finite under h > and can 
be estimated using eq.(2) combined with an appropri- 
ate growth law for the SG short-range order. For a 
sufficiently small h, this teq becomes of an astronomical 
order, and so the blocking of certain modes, spin clusters 
in the present argument, is inevitable in real FCM mea- 
surements. Practically, however, the FCM cusp is a good 



estimate of Tc, i.e., T* 



Tc, if /i is small and the cool- 



ing rate is sufficiently slow. A detailed scenario of 'glassy 
dynamics' for Heisenberg spin glasses is of importance, 
but it may become more complicated due to the presence 
of two types of degrees of freedom: continuous spin and 
(Ising-like) chirality.^'^'^^ It is also of interest to analyze 
the dynamic transition of structural glasses on the basis 
of the present scenario of 'glassy dynamics'. 
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